Necessity of Superposition of Macroscopically Distinct States 
for Quantum Computational Speedup 

Akira Shimiz\J3 and Akihisa Ukena 
Department of Baste Scienee, University of Tokyo, 3-8-1 Kornaba, Tokyo 153-8902, Japan 

Yuichiro Matsuzaki 

Department of Basic Scienee, University of Tokyo, 3-8-1 Komaba, Tokyo 153-8902, Japan and 
Department of Material, University of Oxford, 0X1 3PH, United Kingdom 
(Dated: September 28, 2010) 

For quantum computation, we investigate the conjecture that superposition of macroscopically 
distinct states is necessary for large quantum speedup. Although this conjecture was supported for a 
circuit-based quantum computer performing Shor's factoring algorithm (A. Ukena and A. Shimizu, 
, Phys. Rev. A 69, 022301 (2004)), it needs to be generalized in order to apply to wide classes of 

■ algorithms and/or other models (such as measurement-based quantum computers). To treat such 

general cases, we first generalize the indices for superposition of macroscopically distinct states. We 
then generalize the conjecture, using the generalized indices, in such a way that it is applicable 
unambiguously to general models if a quantum algorithm achieves exponential speedup. Based on 
this generalized conjecture, we further extend the conjecture to Grover's quantum search algorithm, 
\ whose speedup is large but quadratic. It is shown that this extended conjecture is also correct 

. for Grover's algorithm. Since Grover's algorithm is a representative algorithm for unstructured 

problems, the present result further supports the conjecture. 
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I. INTRODUCTION 



We consider quantum speedup for solving computational problems of size L bits, such as the factoring problem (for 
which L is the size of the number to be factored) and the search problem [L is the size of the solution space). In 
the weU-known quantum algorithms of Shor [J 0] and Grover ^ , such problems are solved using quantum computers 
whose number of qubits C > L Since quantum speedup becomes relevant for large L, such quantum computers 

are many-body quantum systems with a large number of qubits C. Since there are many types (and correspond- 
lO ' ing measures or indices) of entanglement for many-body systems interesting to explore which types of 

. entanglement are relevant to large quantum speedup over classical computations |14l421| . 

' This issue has been studied extensively, particularly on Shor's factoring algorithm P, Q and Grover's quantum 
, search algorithm 0. For example, Parker and Plenio demonstrated that the bipartite entanglement as measured by 
the logarithmic negativity is an intrinsic part of Shor's algorithm [l^ . Shimoni, Shapiro and Biham showed that 
^ highly entangled states are generated in both algorithms [isl . [l6| . Oriis and Latorre studied scaling of entanglement 
in three algorithms including Shor's and Grover's [17j |. 

For general algorithms, a few necessary conditions were derived for computational speedup over classical compu- 
^ . tations. Jozsa and Linden showed that for exponential speedup a state is necessary which cannot be factored into 
■ - ' a direct product of states of at most a constant number of qubits [l^ ■ Vidal showed that a necessary condition for 
exponential speedup is that the amount of the bipartite entanglement between one part and the rest of the qubits 
increases with L [l9|. 

We note that one can get a stricter condition by taking the product of these and other necessary conditions, 
which may be obtained by studying other types of entanglement. Such a stricter condition would lead to deeper 
understanding of quantum computations. Hence, it is important to seek for more conditions which are necessary for 
quantum computational speedup. 

As a possible necessary condition, one of the authors conjectured that superposition of macroscopically distinct 
states would be necessary for quantum computational speedup (Ref. [l^l and Sec. IIII[) . Although 'superposition of 
macroscopically distinct states' was only ambiguously defined until recently, clear definition and the corresponding 
index p {1 < p < 2) for pure states were proposed in Refs. according to which a pure state has superposition 
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of macroscopically distinct states if and only if p = 2. The generalization to mixed states was made in Ref. [22], in 
which p is generalized to an index g (1 < g < 2) ; a mixed state has superposition of macroscopically distinct states 
if and only if q = 2. For pure states, p = 2 implies q = 2 and vice versa [22l |. [For mixed states, p is undefined.] 

The entanglement possessed by superposition of macroscopically distinct states is totally different from bipartite 
entanglement which was studied in the previous works ^1§L flij. [Hence, it was called 'macroscopic entanglement' in 
Refs [lllilllllil. However, we do not use this term in this paper because the same term is used in other senses 
by other authors.] For example, some pure states with p = 2 (such as the GHZ state) have very small bipartite 
entanglement whereas some other states with p = 1 (such as energy eigenstates of many-body chaotic systems) have 
very large bipartite entanglement [13, ES]- Therefore, simultaneous requirement (for quantum computational speedup) 
of the superposition of macroscopically distinct states and large bipartite entanglement is much stronger than the 
requirement of either one of them. 

For a circuit-based quantum computer performing Shor's factoring algorithm [HQ, we obtained results that support 
the above conjecture in the previous paper |21|]. It is interesting to study the correctness of the conjecture in other 
algorithms and/or other models (such as measurement-based quantum computers). To explore such general cases 
unambiguously, however, the conjecture needs to be generalized. For example, quantum states in quantum computers 
are not only inhomogeneous but also dependent on instances (i.e., different for different questions of a given problem). 
Since the indices p and q assumed a family of similar states that are spatially homogeneous, the conjecture (which 
was based on p or q) is not strictly applicable to such a general family of states, in its original form. 

Furthermore, since Shor's algorithm is a representative quantum algorithm for solving structured problems [3], it is 
very interesting to study whether the conjecture is correct in the case of quantum algorithms for solving unstructured 
problems. However, the quantum speedup achieved by Grover's quantum search algorithm [3|, which is a representative 
algorithm for unstructured problems 4], is not exponential but quadratic. The conjecture, in its original form, did 
not assume speedup of such a degree. 

In this paper, we first generalize the indices p and q to treat general algorithms and models. We then generalize 
the conjecture, using the generalized indices, in such a way that it is applicable unambiguously to general models 
if a quantum algorithm achieves exponential speedup. Based on this generalized conjecture, we further extend the 
conjecture to the quadratic speedup of Grover's quantum search algorithm. It is shown that this extended conjecture 
is also correct for Grover's algorithm. To show details of evolution of superposition of macroscopically distinct states, 
we also perform numerical simulations of a quantum computer that performs Grover's quantum search algorithm. 

This paper is organized as follows. In Sec.|lTl we generalize the indices p and q. Section [IIII is devoted to generalize 
and further extend the conjecture. Analytic results for Grover's quantum search algorithm are given in Sec. IIVI 
where we will prove that the extended conjecture is correct for Grover's algorithm. We present results of numerical 
simulations of a quantum computer that performs Grover's quantum search algorithm in Sec. |Vl Discussions and 
summary are given in Sec. IVIl 

II. INDICES OF SUPERPOSITION OF MACROSCOPICALLY DISTINCT STATES 

The indices of superposition of macroscopically distinct states were proposed and studied for pure states in Refs. [ill , 
[l^. [23^. and for mixed states in Ref. 22]. To study these indices for states in quantum computers, we here generalize 
their definitions, because, as will be illustrated explicitly in Sees. IIVI and |Vl quantum states in quantum computers 
are not only inhomogeneous but also dependent on instances. Here, an instance is a particular question of a given 
problem. The physical meanings and implications of the indices will also be described briefiy in this section. 

A. Index for a family of pure states 

Let ]7/'^(L))'s be pure states of a system of size L, which are labeled by an index e.g., as 
\ipi{L)),\ip2{L)),\4'3{L)),- ■ ■ . For each value of L, the range of v is given, for example as = 1,2,..., 2^. In a 
quantum computer which solves a decision problem, L corresponds to the size of a certain register, which is usually 
proportional to the size of the input of the problem, and v labels various inputs (see later sections). We do not assume 
that {'tjj^{L)\^p^'{L)) = for ^ . We consider a family F of \^p^{L))'s, 

^'^ {IV.(L)) 1 ah z^'s, all L}, (1) 

which we abbreviate to {\ipviL))}L^^. 

In general, a quantum computer consists of a large number of small quantum systems, such as qubits, which are 
distributed spatially. We call each small quantum system a site, and an operator acting on a single site a local operator. 
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To avoid mathematical complexities, we limit ourselves to the case where every local operator is bounded (i.e., its 
norm is finite). 

Let a(l) be a local operator on site I. We normalize it as ||a(Z)|| ~ 1. Although this normalization condition might 
look too restrictive at first sight, it actually imposes only a weak restriction that makes the maximization operation 
in Eq. ([5]) well-defined [24| , as discussed in Appendix \X\ 

Note that we can use either ||a||£; or \\a\\H as the operator norm ||a||, where 

\\a\\E = sup ||a|V^)|l, (2) 

111^)11=1 



\\a\\H = \^ia,a), (3) 

and (a, b) = Tr(d'''6) denotes the inner product of operators a and b. In fact, both definitions give the same value of 
the index p, because ||a(/)||£; < ||o(/)||// < ■\/d||a(^)||_E and \\a{l)\\H/Vd < ||a(Z)||£; < ||a(/)||fl', where d is the dimension 
of the Hilbert space of a single site. 

By the same symbol a{l) we also denote a{l) which is an operator on the Hilbert space H of the total 

I'm 

system, where !(/') is the identity operator acting on site I' . Using this notation, we define an additive operator A{L) 
as the sum of local operators [ll|, ; 

A{L)=J2&il) (ll«(OII=l)- (4) 
1=1 

Here, we do not assume that a{V) (l' ^ I) is a spatial translation of a{l). 

To simplify notation, we express the expectation value in \'ip^{L)) as We also use the symbols 0,il, to 

describe asymptotic behaviors according to Ref. f3], as summarized in Appendix IBI Furthermore, as described in 
Appendix [Bl a family {/iy(-^)}L,i/ of non-negative functions of L is said to be <d{g[L)) if fy{L) is Q{g{L)) for almost 
every v, i.e., apart from possible exceptional i^'s whose measure (i.e., the number of such z^'s divided by the total 
number of z/'s) vanishes as L goes to infinity. 

For each state consider fluctuation {/S.A{L)^ /S.A{L)) lu , where /S.A{L) = A{L) — {A{L))lu- Its magnitude 

depends on A{L), i.e., on the choice of a(Z)'s. Since ||A(i)|| is upper-bounded, there exists the maximum value, 
max^^^-j (AA(i)^AA(L))Li,. The maximum value is taken for some additive operator A^sx{L), which we call the 
most fluctuating additive operator. Using the maximum value, which depends on L and v, we define the index p of 
the family F (= {\^y{L))} as the positive number (if it exists) that satisfies 

max (V',(i)|Ai(i)tAi(L)|^,(L)) = Q{L-p). (5) 

A(L) 

Note that p does not necessarily exist for a general family. If p exists for a given family, we can show (see Appendix 
El) that 

1 < P < 2. (6) 

When a family {\'(piy{L))} l.v has p — 2 (or p = 1, etc), we also say simply that 'almost every state \tpi,{L)) has p = 2 
(or p = 1, etc).' 

The present definition of p contains those of the previous works, Refs. [lll - fisl . [23j . as special cases. In fact, 
Refs. [illlil,!!^ treated homogeneous states, for which \ip{L+ 1)) was simply defined as the enlarged state of \ip{L)). 
For example, when \^{L)) = |00 • • • 0) -h |11 • • • 1) with L sites, then |V'(i + 1)) = |000 - --O) + ••!) with {L + 1) 
sites. According to the present general definition of p, such a case corresponds to a family whose members for each L 
are identical, i.e., \tp^{L)) = \ip{L)) for all z/. Moreover, Ref. [l3] treated energy eigenstates of homogeneous chaotic 
systems. Since each energy eigenstate is inhomogeneous spatially, this case is different from the case of homogeneous 
states. According to the present definition of p, it corresponds to a family composed of all energy eigenstates in a 
certain energy interval. Therefore, p defined here is a natural generalization of p of the previous works, Refs. [ill - 

The physical meaning of the present p is basically the same as the p of these previous works. That is, almost 
every state of a faniily wi th p — 2 contains superposition of states with macroscopically distinct values of some 
additive operator(s) |lll4l3l . l23l |. We call such an additive operator(s) a macroscopically-fluctuating additive operator(s) 
of the state (or family). For example, when |?/;^(L)) = |00 • • • 0) -f |11 • • • 1) for all v then its macroscopically- 
fiuctuating additive operator is '^i^^zi^), which corresponds to the z component of the total magnetization of a 
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magnetic substance. Since additive operators are macroscopic dynamical variables [ll|, ]23,, two (or more) states 
are macroscopically distinct from each other if they have macroscopically distinct values of an additive operator. 
Therefore, one can state definitely that a state with p — 2 contains superposition of macroscopically distinct states. 

To illustrate how p is useful and necessary for identifying such superposition, we give a few simple examples in 
Appendix [D] 

In the present general definition of p, a macroscopically-fiuctuating additive operator (s) of a given family can be 
different for different L and i^, unlike the case of homogeneous states treated in Refs. [Hi, [H, H^. Therefore, we can 
say the following: For a given family of pure states, if there exists a family of additive operators 

{A^iL) I alliy's, all L} (7) 

such that 

(Ai,(L)tAi,(L))i, = 6(2.2), 

then the family has p = 2. This means that almost all states of the family contain superposition of states that have 
macroscopically distinct values of some additive operator(s). 

It is worth mentioning that there often exist two or more families of such operators, even for a family of homogeneous 
states p^l23j. 



B. Index for a family of mixed states 

The index p, which is defined only for pure states, is sufficient for the concrete analyses in Secs.lIVlandlVl However, 
to state our conjecture in a general form, we need a more general index which is applicable to mixed states. For 
example, measurement-based quantum computers [1^, [2^ often have a larger number of qubits than corresponding 
circuit-based computers. Hence, in such computers, superposition of macroscopically distinct states might appear as 
a mixed state of a subset of qubits \2f\ . Fortunately, generalization of p to mixed states has been made in Ref . [1^ , 
in which a generalized index q was proposed for homogeneous mixed states. We here generalize it to families of more 
general mixed states in order to apply q to wide classes of quantum computers. 

For an additive operator A(L), as given by Eq. and a projection operator f} on TL, satisfying 77^ = 77, we define 
the following Hermitian operator; 

C^. EE [A, [A, f,]] = A^ - 2Af,A + fiA\ (9) 

For a family {pv{L)}l,v of mixed states, which arc not necessarily homogeneous spatially, we define the index q as 
the positive number (if it exists) that satisfies 



L, max(C)L. \=Q{L'i), (10) 

A-.f, J 

where max^ - is taken over all possible choices of A and f). If this q exists for a given family we can show (by slightly 
generalizing the proof in Ref. (22i] ) that 

l<q<2. (11) 

We say a family {pi,{L)} l of mixed states contains superposition of macroscopically distinct states if and only if 
q exists and q — 2. We also say that almost every Pu{L) of such a family contains superposition of macroscopically 
distinct states. For pure states, this is consistent with the corresponding statement based on p, because we can show 
(following the proof in Ref. [23]) that for pure states p = 2 implies q = 2 and vice versa. 

The case of Ref. corresponds to the special case where p^{Lys are homogeneous and independent of 1^ and 
p,j{L + 1) is simply the enlarged state of p^{L). In such a case, q defined here reduces to q of Ref. 2^ . 



C. Properties of states with p — 2 or q = 2 

As application of general theory of Ref. [11], the index p was studied for many-magnon states in Ref. [l3l, for energy 
eigenstates of many-body chaotic systems in Ref. [l3], and for typical many-body states in Ref. J23']. Comparison 
with a measure of bipartite entanglement were also made in these references. Most importantly, many states were 
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found (such as energy eigenstates of a chaotic system) such that they are almost maximally entangled in the bipartite 
measure but their p is minimum, p = 1. Many other states (such as the GHZ state) are also found such that their p is 
maximum, p — 2, but their bipartite entanglement is small. Therefore, the aspect of entanglement which is detected 
by p or g is completely different from that detected by the bipartite measure. Furthermore, it has been shown in 
Ref. [22] that a family of states with q — 2 has a strong L-point correlation, which is Q{L) times larger than that of 
any separable states. Note that any measure of bipartite entanglement cannoi detect such a strong L-point correlation 
for mixed states. On the other hand, the index q does not detect entanglement generated by a small number of Bell 
pairs, whereas measures of bipartite entanglement do. These facts demonstrate that the index q is complementary to 
the measures of bipartite entanglement. 

Note that p is calculated from two-point correlations because fluctuation of an additive operator is the sum of 
two-point correlations; 

L L 

(Ai(L)tAi(i))i, = ^ ^(Aa(O^Aa(Z'))L.. (12) 

1=1 i'=i 

However, this does not mean that p is related only to two-point correlations, because, as mentioned above, pure states 
with p = 2 have q = 2, which means a strong i-point correlation. That is, given the knowledge that a family consists 
of pure states one can say that if the family has p = 2 then it has a strong L-point correlation. 

It was shown in Ref. that p is directly related to fundamental stabilities of many-body states against decohercnce 
and local measurements. Regarding decoherence by weak noises, it was shown that for any state with p = 1 its 
decoherence rate T by any noises never exceeds 0{L). For a state with p = 2, on the other hand, it is possible in 
principle to construct a noise or environment that makes F of the state 0(L^). However, this does not necessarily 
mean that such a fatal noise or environment does exist in real physical systems; it depends on physical situations 
A more fundamental stability is the stability against local measurements, which was proposed and defined 
in Refs. [ll|) HI]- From the theorem proved there, we can say that a state with p = 2 is unstable a gain st local 
measurements, i.e., there exists a local observable by measurement of which the state changes drastically |29| . 

Furthermore, a quite singular property was proved rigorously in Ref. [Tlj : any pure state with p = 2 in a finite 
system of size L does not approach a pure state in an infinite system as L — ^ oo. For readers who are not familiar 
with quantum theory of infinite systems [30| , we give a brief explanation for this fact in Appendix |E1 

These observations indicate that states with p = 2 or q = 2 are quite anomalous many-body states. This led to the 
conjecture of Ref. [lO], which will be generalized in Sec. |in]of the present paper. 

D. Efficient method of identifying superposition of macroscopically distinct states 

Evaluation of a measure or index of entanglement often becomes intractable for large L. Fortunately, this is not the 
case for p because there is an efficient method of calculating p [l^, [l^, [2^ . Since this method assumed homogeneous 
states, we here generalize it to study general families of pure states. 

In this subsection, the dimension d of the local Hilbert space is arbitrary, and we employ \\d\\H defined by Eq. ([3]) 
as the operator norm ||a||. 

Let {bo{l), bi{l), • • • , boil)} be a complete orthonormal basis set of operators on site I, where 6o(0 — l{l)/Vd and 
D = d^ — 1. We expand d{l) as 

D D 

0(0 = 5I^'"^"W' II|cia|' = l, (13) 

where the latter equality comes from |la(/)||// = 1. Let Aa(/) = d{l) — (a(Z)) and Aba{l) = ba{L) — (ba{L)) lv Note 
that A&Q,(Z)'s, unlike &q(/)'s, are not necessarily orthogonal to each other. 

Since A6o(0 = 0: terms with a = do not contribute to AA(L) = Aa(/), i.e., it can be expanded as 

I 

L D D 

Ai(L)=^^Ci„AS„(Z), 5]|qJ2<i. (14) 

1 = 1 a=l a=l 

As a result, X]q=i IczqP 1^ L. Since this normalization is not convenient, we temporarily consider another operator 
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A^' whose expansion coefficients are normalized] 

Ai'(L) = X]EcLAS„(0, X^f]|c;j2^L. (15) 
;=i Q=i 1=1 Q=i 

Here, we do not require that Wia\^ ~ 1 ^^^r every I. The fluctuation of such an operator is calculated as 

L D L D 

(Ai'(L)tAi'(L))i. = EE E E Cc;,,,l^„,p„-(L,^) (16) 

(=1 Q=l !' = 1 a' = l 

for each \tp[y{L)). Here, for a, a' — 1,2, ■ ■ ■ ,D and = 1, 2, • • • , i, we have defined 

Vio^,i,o.'{L,,y) = {Mi{l)^K.{l'))L., (17) 

which can be regarded as elements of a DL x DL Hermitian matrix, which we call the variance- covariance matrix 
(VCM). Therefore, for each \ij^{L)), 

max(Ai'(L)^Ai'(L))i, = e„,ax(i, i^)L, (18) 

A'{L) 

where e^i^^lL, v) is the maximum eigenvalue of the VCM. This and Eq. ([S]) suggest that the following index p^. should 
be useful: 

emax(i,^^) = e(iP=-l). (19) 

In fact, we can show (see Appendix IF|) that ij Pe = 2 then p — 2 and vice versa. For a family of homogeneous states, 
in particular, we can show that p = pe for every value of p (see the last paragraph of Appendix IF|) . Hence, one can 
identify states with p ~ 2 by calculating eniax(ij t^)- This can be done in time Poly(i) because the VCM is a DL x DL 
matrix. 

Furthermore, if Pe — 2 we can find a macroscopically-fiuctuating additive operator from the eigenvector(s) 
Wiamaxi^)} ^hc VCM Corresponding to eniaxiL,^), as follows. If we normalize {c;^„,ax(-^)} 

L D 

EEKmax(i)P=i, (20) 
1=1 a=l 

then the operator 

L D 

Ai-,JL) ^J2T. c;:.nax(i)A6„(0 (21) 

1=1 a=l 

takes the form of Eq. (|15p , and it fiuctuates macroscopically; 

{AAZUl)AAZ,^{L))i^^ = e^UL, i^)L = Q{L'). (22) 
As shown in Appendix |F1 if we put 

C,{L) ^ max [y.^ \cZ.^aAL)\^ , (23) 
then C^{L) — ^{L'^) if pe = 2. Therefore, the following operator (which clearly takes the form of Eq. ([T4|) ) 

L D 



AA^^.AL) cr„max(i)A6a(0, cL„_(L) ^ ^ / , (24) 

V ^}y(L) 



7a max V / 

-lamaxV-^y^^a Hama^\^J = ~ 

1=1 a=l 

also fluctuates macroscopically; 

{AA±,{L)AA'^^JL))l, = e„,,,(i, ,^)L/CAL) = Q{L^). (25) 

We can then construct an additive operator A'^^^^{L) easily from AA'^^^{L), by going from Eq. back to Eq. (|13l) . 
Although A5^a^^(L) is not uniquely determined from AA'^^y.{L) (as discussed in Appendix E)) . this non- uniqueness 
does not cause any difficulty because p is defined only through the fluctuation. 
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III. CONJECTURE ON QUANTUM COMPUTATION 

The conjecture of Ref. [20| is roughly that superposition of macroscopically distinct states should be necessary for 
large quantum speedup. We now generalize it to treat wide classes of algorithms and models. 

We consider decision problems because most computational problems can be reduced, with polynomial overheads, 
to some decision problems. The number of bits or qubits in a computer is allowed to be Poly(L), where L denotes the 
size of the input measured in bits. To be definite, we assume that a quantum computer is composed of qubits (i.e., 
two-level quantum systems) which are separated spatially to each other. That is, if we use the terms of the general 
discussions of the previous section, each qubit is located on its own site. 

We consider the time complexity of problems, allowing both quantum and classical algorithms to have a bounded 
probability of error. In doing so, we assume a classical probabilistic Turing machine as a counterpart of a quantum 
computer. 

A. Exponential Speedup 

To establish notation and to exclude possible ambiguity, we first define exponential speedup, according to convention, 
as follows. 

We consider problems which are not in BPP (Bounded-error Probabilistic Polynomial time) . For a given (decision) 
problem, we denote an instance (input) by i(L, v), where v labels different instances (inputs) of size L. The computa- 
tional time depends not only on L but also on v. For a given classical computer C and for a given quantum computer 
Q, let Tc{L, v) and Tq{L, v), respectively, be computational time for an instance i{L, v), with a bounded probability 
of error being allowed. 

Since we are considering a decision problem which is not in BPP, for any classical computer C there exists a set of 
infinitely many instances that cannot be solved in polynomial time. That is, 

Tc{L,i>) > Poly(L) for some infinitely many instances. (26) 

Here, Tc{L,v) > Poly(L) means that Tc{L^v) — Q,{P{L)) for any polynomial P{L). For a quantum computer Q 
solving such a problem, we say Q achieves exponential speedup if 

Tq{L, v) = Poly(i) for all v. (27) 

According to this definition, Shor's algorithm [l| achieves exponential speedup (if the factoring is not in BPP). On 
the other hand, the Deutsch-Jozsa algorithm fs^] does not achieve exponential speedup because the Deutsch's problem 
is in BPP. 

B. Extra qubits and redefinition of local sites 

In a quantum computer there can exist many qubits (such as ancilla qubits) which are not directly related to 
quantum computational speedup. For example, one can replace a classical circuit that assists a quantum computer 
with a quantum circuit. Then, the size of the quantum computer becomes larger than the original one. It is clear 
that in the enlarged quantum computer only the original part is relevant to quantum computational speedup. 

Therefore, we allow looking only at a subsystem of a quantum computer in order to find out a relevant part, whose 
state (according to our conjecture) would have q — 2. 

Furthermore, one can add extra qubits and circuits to a quantum computer without increasing Tq{L, v) more than 
a Poly(L) factor. For example, to implement quantum error correction fs^, |3^ one can replace each qubit with a 
logical qubit, which is composed of n qubits, where, e.g., n = 9 for the Shor code [s^l- In such a case, correlation 
between two qubits is turned into correlation between two logical qubits, i.e., correlation among 2n qubits. As a result, 
it is possible that a state with g = 2 of the original computer is changed into another state with q <2. However, such 
nonessential decrease of q may be recovered by regarding each logical qubit as a 'local site'. 

Generally, in systems which are composed of discrete sites, a local site (which may be, say, a quantum dot) physically 
has a finite spatial dimension. A set of several neighboring sites also has a finite dimension. Hence, the definition of 
a 'local site' is to a great extent arbitrary. It is therefore possible and reasonable to redefine a set of neighboring sites 
as a new single site |29| . 

In quantum computers, there is further arbitrariness because it is possible to swap the states of two distant qubits, 
paying only polynomial overhead. 

From these observations, we allow all possible redefinition of local sites (accordingly, the number of local sites is 
changed) by regarding two or more qubits, however distant they local site. 
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C. Generalized conjecture 

In order to apply to wide classes of algorithms and models, we generalize the conjecture of Ref. [20] as follows. 

For a decision problem which is not in BPP, consider a quantum computer solving it. If the quantum computer 
achieves exponential speedup, then states having q — 2, whose size is fl{L), appear during computation, for some set 
H of infinitely many instances; 

H = {i{L, I') I some infinitely many {L, i^)'s} (28) 

if 'local sites ' of the quantum computer are appropriately defined. 

This generalized conjecture can be rephrased as follows. After defining local sites appropriately, look at a certain 
subsystem which is composed offl(L) local sites of the quantum computer. Let pk{L,i') be the reduced density operator 
of such a subsystem at the k-th step of the quantum computation. Take some function k^,{L, v), which takes positive 
integral values, of L and v. For some set H of infinitely many instances (Eq. \2^) ). consider the following family of 
states; 

Fk, [H] = {pk, [L, v) I all V and L such that i{L, v) € H}. (29) 

// the quantum computer achieves exponential speedup, one can find an appropriate definition of local sites, a function 
kst{L, v), and a set H , such that q — 2 for the family F]^^ (H). 

We will explain physical meaning of the set H in the next subsection. 

D. Physical meaning of the set H 

If the above conjecture is correct, we can show that H contains some infinitely many instances which cannot be 
solved in polynomial time by any classical computer. That is, H contains infinitely many 'hard' instances. This can 
be seen using reduction to absurdity as follows. 

Suppose that the conjecture is correct but H did not contain infinitely many instances which satisfy inequality 
(PS)) . Then, there would exist a classical computer that solves all instances in H in time Poly(L). By attaching this 
classical computer to Q as a pre-processor, one could obtain another fast quantum computer Q'. However, states 
having q = 2 would not appear in Q' at all, in contradiction to the conjecture. Therefore, H must contain infinitely 
many instances which satisfy inequality (|26p. if our conjecture is correct. 

Note that H is not uniquely determined for a given problem because an appropriate subset of H can be another 
H . To confirm the above conjecture, it is sufficient to find one of many possible iJ's. 

The set H is closely related to a 'complexity core' [s^ - ISTj . A complexity core (or, polynomial complexity core) C 
was defined by Lynch [3^ as an infinite collection of instances such that every algorithm solving the problem using a 
deterministic Turing machine needs more than polynomial time almost everywhere on C. [Note that C is not uniquely 
determined for a given problem because an appropriate subset of C is also a complexity core [s^ - fSTj .] His idea has 
been generalized to complexity classes other than P in Refs. [35- 37]. Intuitively, a complexity core is a set of 'hard' 
instances. The above-mentioned fact shows that H includes a complexity core as a subset. Hence, our conjecture 
claims roughly that states with q = 2 appear for infinitely many instances in a complex core. 

E. Remarks on the conjecture 

Before going further, we make a few remarks. 

The conjecture does not claim that states with q — 2 would be sufficient for exponential speedup; it rather claims 
that they are necessary. Hence, if states with q = 2 appear in some quantum algorithm, it does not necessarily mean 
that the algorithm achieves exponential speedup. 

Moreover, even when a quantum computer does achieve exponential speedup, the conjecture does not claim that 
all states with q = 2 appearing in the computer would be relevant to exponential speedup. In fact, we already showed 
in Ref. [2l[ that although the final state |'0dft) of the computation (before the final measurement) has p = 2 (hence 
g = 2) it is irrelevant to exponential speedup. 

Furthermore, for some states with q — 2 oi size L (such as the GHZ state), one can construct a quantum circuit 
that converts a product state into such a state only in Q{L) steps if the target state with q = 2 is known beforehand 
(i.e., when one constructs the circuit). This fact has nothing to do with our conjecture. The point is that in quantum 
computation the state that appears at, e.g., the middle point of the computation is unknown when the circuit is 
constructed, because the state varies considerably according to the instances. 
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Finally, we note that the results of Ref. [2l| can be understood more clearly according to the present generalized 
conjecture. For example, it was shown in Ref. [2lj that states with p — 2 does not appear when r = 2", where 
n = 1, 2, 3, • • • and r is the least positive integer which satisfies x'' = Ifmod N). [N is a positive integer to be factored, 
and X a random integer co-prime to N which satisfies < x < N [J [2, Hlj.] This fact does not conflict with the 
conjecture. In fact, according to the generalized conjecture, one can exclude the instances with r — 2" from H, or, 
even if one includes such instances in H, such rare instances do not affect the generalized index p of Sec. Ill Al because 
it is defined not by 9 but by 8. 



F. Further extention to Grover's quantum search algorithm 

Shor's factoring algorithm [l| is a representative algorithm for structured problems [4j. For this algorithm, Ref. [2l| 
supported the above conjecture. 

On the other hand, a representative algorithm for unstructured problems is Grover's quantum search algorithm 
Hence, it is tempting to examine the conjecture in Grover's algorithm. However, we cannot apply the conjecture 
(even in the above generalized form) directly to Grover's algorithm because Eq. ([77)) is not satisfied, i.e., it does not 
achieve exponential speedup. Nevertheless, it is often argued that the quadratic speedup of Grover's algorithm is 
significant [?|. Furthermore, Grover's algorithm is known to be optimal, i.e., no quantum algorithm is faster than 
Grover's algorithm by more than a Poly(L) factor in solving the search problem Q. It is therefore very interesting to 
examine the conjecture, if possible, for Grover's algorithm. To make it possible, we here extend the conjecture further 
to Grover's algorithm. Its correctness for Grover's algorithm will be proved in the next section. 

Grover's search problem is the problem of finding a solution to the equation fhix) = 1 among N = 2^ possibilities, 
where fhix) is a function, : {0, 1}^ ^-^ {0, 1}. Let M be the number of solutions, and xi, • • • , xm be the solutions. 
For each L, the solutions specify an instance. That is, (xi, • • • , xm) corresponds to which labels instances as i{L, v). 
According to convention, we regard the number of oracle calls as the computational time. 

In extending the conjecture, we note that the degree of quantum speedup depends on the magnitude of M [see 
Sec. IIVB] . We here consider the case where 

M^0{2"'^) (0<m<l), (30) 

where to is a constant, independent of L. In this case, the degree of quantum speedup is similar to the case of M = 1 
[see Sec. IIVB] . and hence seems significant. 

We extend the conjecture of Sec. IIII CI simply by replacing the condition 'if the quantum computer achieves ex- 
ponential speedup ' with the relaxed condition 'if the quantum computer achieves exponential speedup or it solves the 
quantum search problem using Grover's algorithm for the case of Eq. i30\) \ 



IV. ANALYTIC RESULTS FOR GROVER'S QUANTUM SEARCH ALGORITHM 

In this section, we show that the extended conjecture of Sec. IIII Fl is correct for Grover's algorithm. 



A. Notation 



We first introduce notation. We assume that 2^) ^ M (the number of solutions), because otherwise classical 
computers could solve the problem quickly. 

It seems obvious that an index register, composed of L qubits, is relevant to Grover's algorithm. We therefore look 
only at the index register, although additional quantum circuits would be present in real quantum computers (see 
discussions in Sec. IIII C|) . For each instance i(i, v), where v — (xi, • • • , xm) (see Sec. IIIIFp . we put 

\a{L,v)) = ^ ^ y \x), (31) 

X (/xi,--- ,xm) 

^ 4f E 1^)- (32) 



X { = xi,--- ,Xm 
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Then the state \tpo{L)) just after the first Hadamard transformation (HT) (see Sec. IV A|) is represented as 

6 
= cos - \a{L,v)) +sm - \l3{L,v)), 

where | — >) = (|0) + |l))/\/2, and the angle 9 is given by 



cos- 



N 



Let 0{L, v) be the oracle operator; 



The Grover iteration 



0{L,y)\x) 



-\x) {x^xi,--- ,xm), 
\x) (otherwise). 



GiL,iy)= 2\ML)){ML)\-I{L) 0{L,v) 



(33) 



(34) 



(35) 



(36) 



performs the rotation by angle 9 in the direction \a{L,v)) \l3{L,v)) in the two-dimensional subspace spanned by 
\a{L, v)) and \P{L, v)). The state |?Afe(i, v)) after k (= 0, 1, 2, • • • ) iterations is therefore given by [1] 



\ML,^)) - G{L,vf\i^r>{L)) 
'2k + l 



cos 



|a(L, v)) + sin 



Hence, by repeating the Grover iteration 

R{L) = ( arccos 



2fc + 1 



vr IN 



(37) 



(38) 



times, the state evolves into |V'fl(L)(i, i^)} — |/3(-^, '^)}- Here, \a] denotes the smallest integer among those larger than 
or equal to a. By observing this state in the computational basis, one can find a solution to the search problem with 
probability > 1/2. The range of k is thus 



0<k< R. 



(39) 



B. M dependence of degree of speedup 

According to convention, we regard the number of oracle calls as the computational time. Then, apart from Poly(L) 
factors, 

Tq{L, y) = R{L) = e (v/a7a^) (40) 

for all instances. In contrast, for classical computers there exist infinitely many instances (as L — > oo) such that 

TciL,iy) = e{N/M). (41) 

The quadratic speedup of Tq{L, v) over Tc{L, v) is significant when M is small enough (such as M — Poly(L)). On 
the other hand, no quantum speedup is achieved when M is too large such as M = A^/Poly(L) because then classical 
computers can solve the problem efficiently. 

To be specific, we limit ourselves to the case of Eq. (PO)) . i.e., where M — 0(2™^) [m is independent of L and 
< m < 1], because almost all interesting applications of Grover's algorithm would belong to this case. For example, 
this includes the case of M = Poly(L), whereas the uninteresting case of M = iV/Poly(L) is excluded. 

When M = 0(2"^), we find 

Tc(L,i^) = e (2(1-™)^) (42) 

for some infinitely many instances, and 

TQ(i,i.)-e(V2(i-™)i) (43) 
for all instances. This quadratic speedup seems to be as significant as that of M = 1. 
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C. Family of states to evaluate q or p 



As discussed above, the degree of speedup depends on how the number of solutions M behaves asymptotically as 
a function of L. For clarity, we treat different asymptotic forms of M separately when investigating our extended 
conjecture. 

Suppose that we are given a functional form of A/, which asymptotically satisfies inequality pop . such as M = . 
Then, it will turn out that the set of all instances is an appropriate choice of H of Sec. IIII CI 

H = {i{L,v) I aU L, all v}. (44) 

To construct a family of states Fk, {H), Eq. ([29|) . we specify the number fc* as follows. Since Grover's algorithm simply 
repeats the Grover iteration R{L) times, it seems natural to take /c* — \R{L)/2] . More generally, it seems reasonable 
to take 

h^\RiL)/s-], (45) 

where s is a positive constant independent of L. It will turn out that this choice of is indeed appropriate. A family 
of states Ffc, (H) is thus constructed, for a given fc* , as 

Fk, (H) = {\iPk, (L, ly)) I aU L, all ly}. (46) 

Since all states of this family are pure states, we will evaluate the index p rather than q. 



D. When M = 1 



When M = 1, Eq. dST]) reduces to \a{L,iy)) = IM^)} + 0{1/VN) = \ML)) + 0(1/2-^/2), whereas Eq. 1^ 
reduces to \f3{L, v)) — \xi). Since |V'o(-^)) and |xi) are product states, we find that p = 1 for (the families composed, 
respectively, of) \il}[){L)), \ct{L^ v)) and |/3(L, v)). Hence, p = 1 for (the families of) the initial and final states. 

For intermediate states |V'fe(-^j of interest, it is convenient to investigate Mx = X); ^x{l-) (which corresponds to 
the a; component of the total magnetization of magnetic substances). Note here that in order to show that p = 2 it is 
sufhcient to find one additive observable (which in this case is Mx) which fluctuates macroscopically. From Eq. (|37|) . 
we find 

{i^k{L,v)\Mx\ML,v)) = cosM ^^j L + 0(1), (47) 

{ML,v)\Ml\ML,v)) = cos^ l^-^e^L^ + 0{L). (48) 

Hence, 

{ML,i^)\{^M.?\ML,v)) = isin2((2/c + l)0)L2+O(L). (49) 
For all states in the family F^^ (H) of Eq. (P^ . wc thus find 

(Vfc. {L, iy)\{AMxf\i;k, (L, ly)) - ^ sin^ (^) + 0{L). (50) 
Since s is independent of L, the right-hand side is Q{L^), and thus p = 2 for this family. 



E. When M = Poly(L) 

We now consider the case where M > 2. In this case, unlike the case of M = 1, \f3{L,i^)) has p = 2 for some 
instances. 

For example, suppose that M — 2 and the two solutions for some instance ly are xq = 1010---10 and 
XI = 0101 •••01. Then, p = 2 for \P{L,iy)) = (|a;o) + \xi))/V2 (i.e., for the family {\P{L,iy)) | all L}), because 
{l3{L,}y)\{AMf)^\l3{L,iy)) = n{L^). Here, Mf = J2ii-'^y^zil), which corresponds to the z component of the stag- 
gered magnetization of antiferromagnets. On the other hand, if a;o = 0000 • • • 00 and xi = 0000 • • • 01 for another 
instance v' then p = 1 for \(3{L,iy')) = (|xo) -I- 1x1))/-^^. 
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Therefore, when M > 2, p of the final state depends on the instance, i.e., on the nature of the solutions. (On the 
other hand, it is clear that p = 1 for the initial state.) 

To compute p of states in intermediate stages of computation for M > 2, we first consider the case where M = 
Poly(L). In this case, Eq. still holds for {ipkiLji')) with fc > 1, and thus the discussion following Eq. also 
holds. Therefore, we again find that p = 2 for the family Fk, {H) of Eq. p6|) . 

We can obtain the same conclusion when M = 0(2^ ), where k is a constant independent of L and < k < 1. 
Instead of showing this, we shall derive the same conclusion when M is even larger in the next subsection. 



F. When M = 6(2™^) 

We now study the case where 

M = 6(2™-^) (0 < TO < 1), (51) 

where to is independent of L. This is the upper limit of M that satisfies condition ([5(7|. 

As mentioned above, p of the final state depends on the nature of the solutions (whereas p = 1 for the initial state) . 
Since this might not be trivial when M is as large as 8(2™^), we give an example for M — 2^/^ = VTV. Suppose 
that the solutions for some instance v are as follows; Xj — j for j = 1, 2, . . . , whereas xj^s for j > + 1 are 

muhiples of 2%/iV (less than N). Then 



- -L|00...0^^...^) + -i=|^^...^00-..0)+O(l/V]V). (52) 

Since {/3{L,i^)\{AM' ^)^\P{L,iy)) = 0{L^), where M'^ = Efif '^^^(O + Ef=L/2+i i(0> we find that p = 2 for this state. 
By contrast, if xj — j for all j for another instance v' , then |/3(L, i/')) — |00 • • • — • • • -^), for which p = 1. 

To compute p of states in intermediate stages of computation, we note that \a{L,i')) = |?/'o(-^)) + 0(2^^^^™)^/^) 
when M = 9(2™-^). Hence, abbreviating l^pkiL,^)), |a(L, z/)), |/3(i,i')) to IV'fe), \ct), \I3), respectively, we find 

ii^klMM = cos^ (a|Af,|a) + sin^ (^^^0^ (/3|M,|/3) + 0(1), (53) 

which yield 



(V-fclMjlV-fe) - cos^ [-^^ I { ] {mM+0{l), (54) 



(V'fe|(AM,)2|^,) = cos2 f ^^^j (a|(AA4)'|a) +sin2 f (/?|(AM,)2|/3) 

+ isin2((2fc + l)0)((a|M,|a)-(/?|M,|/3))' + O(L). (55) 

As shown in Appendix[Gl {a\Mx\a) = L + 0{1) and {(3\Mx\(3) = KL, where K is independent of L and < K <1. 
Hence, for h of Eq. gSJ, we find {iljk,\{AMx)^\tpk,) = 0{L^ )- Th erefore, p^2 for the family Fk, (H) of Eq. pH]). 
We have thus proved that the extended conjecture of Sec. IHIFl is correct for Grover's algorithm. 



G. Intermediate values of p 



As the quantum computation proceeds (i.e., as k is increased), p increases from 1 for the initial state |?/;o(i)) to 2 
for (L, i^)). In the transient steps, p takes intermediate values between 1 and 2. 

Such intermediate values are taken also, e.g., by states of quantum many-body systems at critical points of continu- 
ous phase transitions, where two-point correlation functions decay according to power laws as a function of the distance 
between two points. Hence, one might expect some universal properties of p, as critical exponents in continuous phase 
transitions have. 
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For states of quantum computers in the transient steps, however, the intermediate values of p are not universal; 
they depend on details such as the nature of the solutions. Since we are not interested in such non-universal values 
of p in this paper, we have focused on the universal result, which is directly related to our conjecture, that p — 2 for 
the family composed of (i, '^))'s. 



V. EVOLUTION OF QUANTUM CORRELATIONS IN GROVER'S QUANTUM SEARCH ALGORITHM 

As summarized in Sec. [Ill the index p is calculated from two-point correlations of local operators. For pure states, 
p = 2 if they have two-point correlations of 8(1) between O(L^) pairs of sites. [As discussed in Sec. Ill C| for pure 
states such strong two-point correlations imply a strong L-point correlation.] As discussed in Sec. IIIDi the existence 
of such correlations can be detected by the asymptotic behavior (as L —> oo) of the maximum eigenvalue Cmax of the 
VCM, because, roughly speaking, Lemax is proportional to the number of pairs of sites whose correlation is of Q(l) 

M- 

Although p is simpler and more convenient for stating the conjecture, Cmax has more detailed information about two- 
point correlations. [For example, a state with Cmax = L has stronger two-point correlations than one with Cmax = L/IO, 
whereas both states have p = 2.] It is therefore interesting to study how emax evolves as the computation proceeds 
from a small value (corresponding to p = 1) at the initial state to larger values. It describes how two-point correlations 
evolves (until a strong L-point correlation develops for p — 2, as discussed in Sec. Ill Cl) . 

In this section, to investigate the evolution of Cmax, we simulate numerically a quantum computer that performs 
Grover's algorithm. In the simulation, we study more states than those studied in the previous section, where we have 
studied \4>k{L, v)) = G^{L, j/)|'0o(L)) (A: = 0, 1, 2, • • • , R{L)). In actual quantum computations, the Grover iteration 
G{L,i') may be realized, e.g., as a series of local and pair-wise operations [4]. Hence, many intermediate states 
appear between the computational steps corresponding to \4'kiL,i>)) and [i/ife+i (L, i^)). Although it was sufficient to 
investigate \ipk{L,vys for confirming the conjecture, we also study such intermediate states to see more details. 

We simulate two cases, M = 1 and M = 2, because these cases are most fundamental. The solution(s) xi (and X2) 
is chosen randomly. We have confirmed that this random choice of a solution(s) makes no significant difference on 
the results of the numerical simulations presented below. 



A. Formulation of simulation 



We explain our simulation for the case of M — 1. Simulation for M — 2 has also been performed similarly. 
Since M — 1, we can simply take — xi. As in the previous section, we consider the index register composed of L 
qubits. The register is initially set to be in the following product state; 

IV'init(i)) - I00---0). (56) 

Firstly, the HT is performed by successive applications of the Hadamard gate on individual qubits, and the quantum 
state evolves into \iPq{L)) of Eq. ([55)1 . Then we apply the Grover iteration G{L,i'), which consists of two HTs, an 
oracle operation 0(L, v), and a conditional phase shift P{L) [1]; 

[-m {x>0). 

Each HT requires L operations of the Hadamard gate. The oracle 0{L, v) requires its own workspace qubits 
and computational time. However, since the oracle is not a proper part of the Grover's algorithm, we simulate the 
operation of 0(i, v) as a one-step operation, and its workspace is not included in the simulation. The execution of 
P{L) requires ^1{L) pairwise unitary operations. For simplicity, however, we simulate P{L) as a one-step operation. 
Hence, each Grover iteration is simulated by 2L -|- 2 steps of operations. 

After application of the Grover iterations R{L) times, the state |V'o(i)) evolves into 

G«(^)(L,z.)|Vo(i)) = \i^RiL){L,v)) ~ \P{L,v)) = (58) 

Finally, by observing this state one can obtain the solution xi with a sufficiently high probability. We do not simulate 
this measurement process. The total computational time (steps) Tq(L, v) in our simulation is thus 



Tq(L, v)^L+{2L + 2)R{L) ^ (d{LV¥) 



(59) 
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for all V. This is larger than Tq{L, v) of Sec. llVl onlv by a polynomial factor. For each instance Tq{L, v) different 
states, including |-0fe(L, i/))'s of Sec. llVi appear during computation. 

In our conjecture, we have allowed (i) looking only at a subsystem and (ii) all possible redefinition of local sites. For 
the present model of a quantum computer, however, we have confirmed (from the following results and the results of 
the previous section) that they are unnecessary. That is, in the present model, we can confirm the extended conjecture 
by simply calculating the index p of states of the index register. 

To find states with p = 2, we calculate eniax(i, v) of Sec. IIIDI By plotting the L dependence of eraax{L, v), we can 
determine the value of Pe through Eq. p^ . When states with = 2 are found, they have p = 2 because, as discussed 
in Sec. Ill Dl ii Pe — 2 then p — 2 (and vice versa). We will also plot how eniax(-^, i^) (for fixed L) grows and decays as 
the quantum computation proceeds because it is instructive and impressive. 

In defining the VCM of Eq. (fTT]). we take ba{l) = <7a{l) (Pauli operator on site I and a — 1,2, 3), i.e, 

Vla.,l,a.'{L,l^) = {Aa^{l)Aaa.,{l'))Lu- (60) 

In the following, for the sake of conciseness, we will often denote eniax(^, i'), Via^i'a'iL, j/), and so on simply by emax, 
Via^i'a' , and so on. 



B. Results of simulation for M — 1 



Figure [U plots evolution of e^ax for L = 8, 9, 10, 12, 14 and M = 1, when xi = 19, 388, 799, 1332, 9875, respectively. 
It is seen that on the whole the curves are exponentially expanded along the horizontal axis as L is increased. 

Figure[2]is a magnification from the 1st to 40th steps for L — 8, whereas a magnification from the 1005th to 1155th 
steps for L = 14 is shown in Fig. [3] 
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FIG. 1: The maximum eigenvalue emax of the VCM of quantum states appearing in Grover's quantum search algorithm for 
L — 8,9, 10, 12, 14 when xi = 19, 388, 799, 1332, 9875, respectively, as functions of the step of the algorithm. The horizontal 
line represents the value of emax for product states, emax ~ 2. 

It is seen that Cmax — 2.00 for all states from l-fAinit) to |V'o)i i-C, during the initial HT (from the 1st to 8th steps 
in Fig. [21 denoted as 'HT'). This is because all these states are product states, for which we can easily show that 
Gmax — 2 f Appendix [HI . When the stage of Grover iterations begins, emax grows gradually, as seen from Figs. [1] and [2] 
In each Grover iteration. Figs. [2] and [3] show that Cmax changes when the oracle operator O is operated, whereas it is 
kept constant during the subsequent HT. Then, it changes again when P is operated, whereas it is kept constant again 
during the subsequent HT. As the Grover iterations are repeated, emax (hence, quantum correlation) continues to 
increase as a whole, until it takes the maximum value after about R{L)/2 times applications of G. Further applications 
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FIG. 2: A magnification of Fig. [T] from tfie 1st to 40tli steps for L = 8. Here, O and P represent the oracle operation and 
conditional phase shift, respectively, in a single Grover iteration G. 
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FIG. 3: A magnification of Fig. [T] from the 1005th to 1155th steps for L — 14. Here, O and P represent the oracle operation 
and conditional phase shift, respectively, in a single Grover iteration G. 



as seen from Fig. Il] toward e„ 



2.00 for \ipB.(,L))y which is approximately a product state as seen 



of G reduce emax, 
from Eq. (|55|) . 

These results indicate that our construction of the family of states, Eqs. (j45|) and (|46l) . is natural. Although we 
have already shown in Sec. |IV| that p = 2 for such a family, it is instructive to plot ema.^{L,v) of |-0fe, (L, i^)) as a 
function of L when the solution Xi is randomly chosen (i.e., an instance u is randomly chosen) for each L. Figure 
IDplots Binax's of \ip\R/2\), |'i/'[fl/3] ) and |'0[fl/4] ) for such randomly-chosen xi as functions of L. [We have confirmed 
that almost identical curves are obtained for other choices of Xi as well.] Since Cmax's tend to be proportional to L 
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for large L, we can confirm that p = pe = 2. 
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FIG. 4: fimax's of lV'ffl/2i) (diamonds) , [V'ffl/si) (crosses) and \ipiR/4,]) (squares), as functions of L. They all show that 
p — Pe = 2. The dotted lines are the guides to the eyes, whereas the horizontal line represents the value of Cmax for product 



C. Results of simulation for M = 2 



When M = 2, the result of Sec. IIV El indicates that p = 2 for most states in the Grover iteration processes whereas 
p of the final state depends on the nature of the solutions. Figure [5] illustrates this clearly, where we have plotted 
evolution of Cmax for two cases, p — 2 (profile 1) and p — 1 (profile 2) for the final state. For \^[r/2]), \fp[R/3]) and 
|V'[_R/4])j on the other hand, we obtain results similar to Fig. |3]for both cases. Hence, p = 2 for these states. This 
result visualizes how our extended conjecture holds when M > 2. 



VI. DISCUSSIONS AND SUMMARY 



We have studied the conjecture that superposition of macroscopically distinct states would be necessary for sig- 
nificant speedup of quantum computers over classical computers. This conjecture was previously supported for a 
circuit-based quantum computer performing Shor's factoring algorithm. To treat general algorithms and models, we 
have generalized the indices p and q for superposition of macroscopically distinct states. We then generalize the 
conjecture in such a way that it is applicable unambiguously to general models if a quantum algorithm achieves an 
exponential speedup. We further extend the conjecture to the speedup achieved by Grover's quantum search algo- 
rithm. This extended conjecture is proved to be correct for Grover's algorithm. Since Grover's and Shor's algorithms 
are representative ones for unstructured and structured problems, respectively, the present results and the results 



of Ref. 21[ support strongly the conjecture. To see details, we have also presented, by numerical simulation, how 



quantum correlation evolves and the superposition of macroscopically distinct states develops as the computation 
proceeds. 

Jozsa and Linden previously showed that entanglement over a cluster whose size is larger than 0(1) is necessary 
for an exponential speedup [18|. For p — 2, on the other hand, entanglement over a cluster whose size is larger than 
0{L) is necessary [111 Il2|. The present conjecture imposes a stronger condition in this sense. 

Moreover, Vidal showed that a necessary condition for an exponential speedup is that the amount of the bipartite 



;peed 

entanglement between one part and the rest of the qubits increases with L [l^ . The bipartite entanglement which was 
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FIG. 5: The maximum eigenvalue emax of tlie VCM of quantum states appearing in Grover's quantum search algorithm for 
L = 10 and M — 2, when xi = 2,X2 = 1023 (profile 1) and when xi = 511, X2 = 512 (profile 2), as functions of the step of the 
algorithm. The horizontal line represents the value of emax for product states, emax = 2. 

studied by him is totally different from the entanglement possessed by superposition of macroscopically distinct states. 
Therefore, the simultaneous requirement (for speedup over classical computations) of Vidal's condition and and the 
present conjecture is much stronger than the requirement of either one of them. That is, for quantum speedups both 
superposition of macroscopically distinct states and a sufficient amount of the bipartite entanglement are necessary. 

It is also interesting to explore relation between our results and the problem of time-optimal quantum evolution 
[39I lioj . In the latter case the optimal evolution requires large fluctuation of Hamiltonian, while in the former fast 
quantum computations require states with large fluctuations of additive operators. This suggests a possible relation 
between the quantum speedup and the optimal evolution j40| . 

Acknowledgments 

We thank M. Koashi and A. Hosoya for valuable discussions. This work was supported by PRESTO, Japan Science 
and Technology Corporation, by Grant-in- Aid for Scientific Research No. 18-11581, and by KAKENHI No. 22540407. 

Appendix A: On the restriction that ||a(OI| ~ 1 

Suppose that A = X)/ '^(O an additive operator. In the present paper, we have required that ||a(OII = 1- To 
understand technical details about this, the following examples would be helpful. 

ex.1 The operator A'l = a{l)/2 is not an additive operator according to the present definition, because the norm 
of the local operators is not unity. However, by simply multiplying A'l by 2, we can obtain an additive operator 
Ai = 2A'l {— A). Fluctuations of Ai and A'l differ only by a constant factor. 

ex.2 The operator A2 = X]!=odd ^(0 is not an additive operator according to the present definition, because the 
norm of the local operators for even I vanishes. However, this operator has the same fluctuation as the additive 
operator A2 = E;=odd "(0 + E/=cvcn UO- 

ex.3 The operator A^ = ^1 + '^"Y^) '^(O '^''^ additive operator according to the present definition. However, 
there always exist real numbers a, /3{l) such that 03(1) = a (^l + '^(0 + /3(0l(0 has unit norm for every 

Then, A3 = J^i ^3(0 additive operator, and fluctuations of A3 and A3 differ only by a constant factor. 



profilel 




profile2 

: ]..- "\ 
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Therefore, operators like A'(,A2,A'^ are essentially included (as Ai, A2, A3, respectively) when taking max^j.^^ in 
Eq. ©. ^ ^ 

The point is that one can modify A", A2, A';^ in such a way that fluctuations of the modified operators Ai,A2, A3 
(which are additive operators) have the same order of magnitude as those of A", A2, A3. Note that this modification 
is noi unique. For example, from A'{ one can also construct an additive operator 2A" — a(l) + 1(1) = 1(1) +X]i>2 ^(Oj 

which has the same order of fluctuation as Ai . This non- uniqueness does not cause any difficulty because p and q are 
defined by the order of magnitude of fluctuations. 

Appendix B: asymptotic notation 

Let f{L) and g{L) be non-negative functions of a positive variable L. Following Ref. we use the following 
asymptotic notation: 

/(i) = 0(.9(L)) ^ /(i) < KgiL), (Bl) 

f{L) = n{g{L)) ^ Jg{L) < f{L), (B2) 

fiL) = e{giL)) ^ Jg{L) < f{L) < Kg{L), (B3) 

as L — > 00, where J, K are some positive constants. 

Let fu{Lys be non-negative functions, which are labeled by an index of a positive variable L. We consider a 
family that consists of these functions, i.e., a family of real values, {/i/(i)}L,iy. Assuming that the number of possible 
values of v increases to inflnity as L — ^ 00, we use the following asymptotic notation: 

f^{L) = 0{g{L)) UiL) = Oig{L)) for almost every z/, (B4) 
U{L) = Ti{g{L)) ^ /^(L) = rj(c/(L)) for almost every I/, (B5) 
U{L) ^ e{g{L)) ^ U{L) = e(5(L)) for almost every i/. (B6) 

For example, if = 1, 2, . . . , 2^ for each L and /^(L) = L'^{1 + sinL)(l/L + l/L"), then f^{L) = Q{L'^), whereas 
it is not 6(^2). 

Appendix C: Range of p 

In this appendix, we show that 1 < p < 2. As the operator norm ||a||, we employ ||a||_E, deflned by Eq. ([2]), in this 
Appendix. Since (^lS.A{L)^ ^ < (^A{L)^^ < L^, we flnd that p <2. To prove that p > 1, we use the following 
Lemma; 

Lemma: For any state, which can be a mixed state, there always exists an local operator a^,(l) which satisfles 

||a*(/)||B = 1 and {Aa^l)^ Aa^l)) = 1, (CI) 
Proof: For a given state p{L), its local density operator pi = Tr;/(-^/-)p(L) can be diagonalized as 

d d 

pi = Y,w,mi){jj\, E^^w^i' (C2) 

j=i i=i 
where {\j, l)}j is a complete orthonormal set of site I. Take 

d,{l)= {\j,l){j + l,l\ + h.c.) (where (d+l,/| = (l,/|), (C3) 

J— odd 

which is a non-Hermitian local operator with |la*(Oll-E — 1- Since Tr [p(L)a^{l)] = and 

d 

a*(0^a*(0=El^'^)O',^l = i(0, (C4) 

we flnd that (Aa* (O^Aa* (?)) = 1. ■ 
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Using this Lemma, we now show the following theorem, from which it is obvious that p > I. 
Theorem: For any state, which can be a mixed state, there always exists an additive operator which satisfies 

(Ai(L)tAi(L)) > L. (C5) 

Proof: We use the induction method. We define 

k 

i(fc)^^a(0 (116(011 = 1), (C6) 
1=1 

ioT 1 < k < L. When k — 1, the above Lemma shows that there exists A{1) such that {AA{1)^ AA{1)) > 1. Now, 
assume that there exists A{k) such that {AA{ky AA{k)) > 1. From the above Lemma, there exists a local operator 
a«(fc + 1) on site k + 1 such that ||a*(fc + 1)|| = 1 and (Aa*(fc + l)'''Aa*(fc + 1)) = 1. So, construct A{k + 1) as 

i(fc + l) =i(fc) + a»(fc + l). (C7) 

Then, 

{AA{k + l^AAik + 1)) = {AA{k)'' AA{k)) + \^{AA{k)'' Aa^k + 1)) + c.c.} + 1. (C8) 

Therefore, if |(Ai(fc)tAa*(fc + 1)) + c.c.j > then (Ai(fc + l)tAi(A; + 1)) > k + 1. If, on the other hand, 
|(A4(fc)tAa*(fc + 1)) + c.c.j < then reconstruct A{k + 1) as 

A{k + l) = A{k)~a4k + l). (C9) 
This gives {AA{k + 1)1' Ai(fc + 1)) > A; + 1. ■ 



Appendix D: Simple examples of states with p = 2 and p = 1 

To explain why p is useful and necessary, we give a few simple examples of homogeneous states. 

For the GHZ or 'cat' state, |V'cat(.^^)) = (|00---0) + |11 • • • l))/-\/2, a heuristic discussion would be sufficient to 
judge that it contains superposition of macroscopically distinct states, because it is simply superposition of two 
macroscopically distinct states and their coefficients do not vanish as L — >■ oo. 

By contrast, heuristic discussions will be ambiguous for the following state; 

l^dw(i)) = -7^= (|100 • • • 00) + |110 • • • 00) + |111 • • • 00) + • • • + |111 • • • 10)) , (Dl) 

where 'DC stands for domain wall (see below). This is superposition of L states, each state differs from the preceeding 

state only in a single qubit, and the weight of each state vanishes in the _L — > limit. 

Even for such a state, by using the index p we can easily judge whether superposition of macroscopically distinct 
states is contained. In fact, we can show that p = 2 for lipdw), hence it does contain superposition of macroscopically 
distinct states. Physically, this state is superposition of single-domain-wall states. That is, let + 1) be the state 
which has a domain wall between sites j and j + 1, then 

1 ^"^ 

\i;d^{L)) = -^J2\j;j + l) (D2) 

Since the domain wall is clearly a macroscopic object, this state contains superposition of macroscopically distinct 
states. This can also be understood by rewriting IV'dw) as 

. L/2-1 L-1 

l^dw(i)) = -y== J2 li;i + i) + ^== E + (d^) 

= ^= [a state whose domain wall is located in the left half] 
v2 

-|-— ^[a state whose domain wall is located in the right half]. (D4) 
v2 
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The two states in this expression are clearly distinct macroscopically. 

On the other hand, the 'W state' \ip^{L)) = (|1Q0 • • • 0) + |010 • • • 0) H h |000 • • • 1)) has p = 1, hence it 

does not contain superposition of macroscopically distinct states [HI, [12] ■ This is reasonable because the W state 
corresponds to normal states in condense-matter physics, such as a Frenkel exciton excited in an insulating solid f2l| . 

Since both the W state and product states have p = 1, they belong to the same class with respect to superposition of 
macroscopically distinct states. In contrast, states with p — 2 are extremely abnormal in view of many-body physics, 
as discussed in Sec. Ill Cl and Ref. 

Regarding the index q for mixed states, see Ref. [l^ for several examples. 



Appendix E: Approach of a pure state with p = 2 to a mixed state as L — !> cxd 

It was proved rigorously that any pure state with p = 2 in a finite system of size L does not approach pure states in 
an infinite system as L — >■ oo pj| . Although this might sound strange to the reader who is not familiar with quantum 
theory of infinite systems [131, its physics can be understood as follows. 

In quantum theory of finite systems, all possible representations are equivalent (unitary equivalence). In quantum 
theory of infinite systems, in contrast, many representations can exist which are not equivalent to each other [soj . 
Among two or more inequivalent representations, we have to choose one that is suitable for describing the physical 
states of interest [s^l- This makes the quantum theory of infinite systems much different from that of finite systems. 
The above strange fact comes from this great difference. 

As the simplest example, consider a cat state \ip{L)) = (|00- • • 0) + |11 • • • 1))/V2 of size L. If L is finite, there 
exist observables that have non- vanishing matrix elements between |00 - • - 0) and |11 • • • 1). The expectation values 
of such observables discriminate between the pure state {ipiL)) and the mixed state p{L) = (|00 • • • 0) (00 • • • 0| + 
|11 • • • 1)(11 • ■ ■ l|)/2. If we take the L ^ oo limit, however, the quantum theory of infinite systems requires that every 
observables should be a function of field operators within a finite region in an infinite space [30]. As a result, there 
are no observables that have non- vanishing matrix elements between the L ^ oo limits of |00 • • • 0) and |11 • • • 1). 
This implies that limL^ao is not a pure state. Here, the rigorous definition of pure states in Ref. [s^l are used 

instead of the (over)simplified definition = p, because the latter can be used only for (an irreducible representation 
for) finite systems. 

More mathematically speaking, if \il){L)) has p = 2 then limL_i.oo |'0(^)) is not a vector state of an irreducible 
representation. For details, see Ref. [ll| and references cited therein. 



Appendix F: Equivalence of pe = 2 and p = 2 

In this appendix, we show that if — 2 then p = 2 and vice versa. For simplicity, we will omit V and '(i)', i.e., 
we will abbreviate A'{L), \ip^{L)), C^{L) and so on to A' , {ip), C and so on, respectively. We first note that 



aA'\ 



L D 



^5^cLAS„(/)|^) 



1=1 a=l 



< 



L D 



EE 

.1=1 a=l 



L D 



5^5] aS„(z)|V') 



.1=1 a=l 



< LOiL), 



(Fl) 



from which p^ < 2. Here, we have used the inequality J^k l^j^jlljk < 12k^j l-'^jf Ivj'^f — 

l^jl^jl^ '^j'12k\yj'k\^ ' which holds for arbitrary complex numbers Xj,yjk. On the other hand, it is clear from 
the definitions that 



P < Pe- 



(F2) 



Therefore, if p = 2 then pe = 2. 

To show the inverse, we assume that pe = 2, i.e., Eq. ([22) holds. Without loss of generality, we also assume that 
niax,_„ |c;„j„a,xl^ = Ic'iimaxl^ • I* is clear from Eq. ^ that 



< L. 



(F3) 



Our purpose is to show that 



Ic'n, 



0(L° 



(F4) 
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because it yields C — 0{L^), which gives Eq. ([25]) (implying p = 2). Equation (|22)) can be rewritten as 
D f L D \ L D 

^ ^ C^Q, niax'^lQ' max^la,la' + I ^ ^ ^ ^ "^Za max'^la' max^a,la' ~^ C.C. j + ^ ^ ^ ^ '^la ma.x^l' a' ma.x^la,l' a' — 0{L ) 
a,a' = l y/=2 q,q' = 1 J l.l'=2a.a'=l 

(F5) 

The first term is positive (because the VCM is a non- negative Hermitian matrix) and < 0{L). The second term is 
estimated as 



L D 



'^la max'^la' max^a,la' 



1=2 a,a' = l 



L D 

< DWii 



Wr'* 

/ J / J ^la max 

1=2 a=\ 



max 



< i?|c'n„,,J V^(£ - 1) VI max \Vio.,iA 

< |c'n„_|0(L)<0(i3/2), (F6) 



where we have used inequality (jF3p and |x • 1| < |x||l| = which holds for an arbitrary n-dimensional vector x. 
Therefore, the third term of Eq. (|F5[) should be 



L D 



E E ^'l*an..Aa'n...Vl^,l'a'^OiL^). (F7) 
l,l'=2 a,a' = l 

On the other hand, let 

max' 

(F8) 

a max 

then the following operator, which does not involve an operator on site 1, takes the form of Eq. (1151) : 

L D 

Ai'^^^cLA6„(0. (F9) 

1=2 a=l 

We therefore have, using Eqs. (|F6P and (|F7|) . 

< (A<LAi-,Ji.-(Ai-tAi-)^, 

_D / L D 

— y ] ''let max'^la' max^la,!"' "'^ I ^ ^ ^ ^ ''ia max'^la' max^a,!"' + CC. 
Q,a' — 1 \ ^— 2 a.a' = l 



11 max I 



|2 



L D 



E E 

^la max '-■/'a' max^o:,/''a' 



L 

i,i'=2 a,a' = l 

< |c'n_,pO(L0) + |c'n_JO(L) - '""^ '' o(L^) 

< Ic'iimaxI [|c'n^,,|0(L°) + 0(L) - |c'ii^,,|0(L)] . (FIO) 

Therefore, Eq. (jF4p should be satisfied because otherwise the last line would become negative in contradiction with 
the first line. 

For a family of homogeneous states, in particular, p — Pe for every value of p, because in this case the VCM has 
the translational invariance and thus an eigenvector corresponding to Cmax is also translational invariant. [Even when 
Gmax is a degenerate eigenvalue, one can construct a translational-invariant eigenvector by taking linear combination of 
eigenvectors corresponding to e,„ax-] From Eq. (pil)) . this means C = 0{L^) for the additive operator that is composed 
of such an eigenvector, and thus p > Pe- From inequality (jF2p . this yields p = Pf.. 



Appendix G: Proof of {/3|M^|^) = fcL (0 < < 1) for M = 2"^ 

In this appendix, we prove (/3|Ma;|/3) = KL {Q < K < \) for M = 2™-^ where K is independent of L. We use the 
following lemma: 
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Lemma: There exists a real number K independent of L such that 2^^^ — KY ^ > 1 and \ < K <1. 

Proof. Consider the foUowing function; 

/(fc) = (l-m)log2 + A;logfc+(l-fc)log(l-A;) < A; < 1), (Gl) 

where < m < 1. Since f{\) = — mlog2 < 0, f{k) — )• (1 — m) log2 > (A; 1), and the function f{k) is continuous 
at i < fc < 1, we can apply a intermediate- value theorem and therefore there exists a real number K such that 
f{K) > and 1/2 < if < 1. ■ 

Now, we prove (^|Mx|/3) = fci (0 < A; < 1) for M = 2™^. Let P(Mx = -L + 2j) be the probability of getting a 
value —L + 2j when one performs the measurement of where j = 0,1, - ■■ ,L. It is represented as 

P{M, = -L + 2j) = J2\{M, = -L + 2j, r < (^) niax[ | (M, = -L + 2j, i.|/3) f ], (G2) 

where \Mx = L — j, v) is an eigenstate of M^, and v = 1,2, - ■ ■ , (^) labels degenerate eigenstates. Since 

1 ^ [m 
I (M, = -L + 2j, v\p)\< y^\{M, = -L + 2j, v\xn) \<\h^, (G3) 



n=l 

we have, using Stirling's formula n! ~ V27rn(-)", 



P(M. = -i + 2i) < f ^) 1^ - , ^ ^ — (G4) 



From the above Lemma, there exists a real number K independent of L such that 2(^~'"^(1 — Ky~^K^ > 1 and 
^ < K < 1. Hence, 

L"P{M^ = -L + 2j) ^0 {KL<j<L) as L oo, (G5) 
where n is independent of L. This yields 

KL-l 

^ E + = -L + 2i) (G6) 

as -L DO, from which we conclude that {j3\Mx\P) < KL for large L. ■ 



Appendix H: Maximum eigenvalue of the VCM for product states 



In this appendix, we show that Cmax = 2 for a product state, \'^) = \4>i)ii where \4>i)i denotes a state of the 

qubit at site / (= 1, 2, • • • , V). The VCM of such a state is block-diagonal; 



/Vi ••• 

F2 



\ 







(HI) 



V Vl / 

where is a 3 x 3 matrix whose a/3 element (a, (3 — x,y, z) is given by 

{ViU - {M^aimpim) - {cf>i\a^{l)\cbi){cj,i\ap{lMi) 



(H2) 



Therefore, Cmax is given by the maximum one among the maximum eigenvalues of Vj's. By a unitary transformation 
of this 3x3 matrix such that becomes an eigenstate of the transformed (t^(Z)j we can transform Vi into 



Vi = 




(H3) 
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Since the maximum eigenvalue of this matrix is 2, we find that Cmax = 2. 



[1] P. W. Shor, in Proceedings of the 35th Annual Symposium on the Foundations of Computer Science, edited by S. Goldwasser 

(IEEE Computer Society, Los Alamitos, CA, 1994) p. 124. 
[2] A. Ekert and R. Jozsa, Rev. Mod. Phys., 68, 733 (1996). 
[3] Lov K. Grover, Phys. Rev. Lett. 79, 325 (1997). 

[4] M. A. Nielsen and L L. Chuang, Quantum Computation and Quantum Information (Cambridge University Press, Cam- 
bridge, 2000). 

[5] C. H. Bennett, S. Popescu, D. Rohrlich, J. A. Smolin, and A. V. Thapliyal, Phys. Rev. A 63, 012307 (2000). 

[6] A. Miyake and M. Wadati, Phys. Rev. A 64, 042317 (2001). 

[7] D. A. Meyer and N. R. WaUach J. Math. Phys. 43(9), 4273 (2002). 

[8] J. K. Stockton, J. M. Geremia, A. C. Doherty, and H. Mabuchi, Phys. Rev. A 67, 022112 (2003). 
[9] O. F. Syljuasen, Phys. Rev. A 68, 060301 (2003). 
[10] F. Verstraete, M. Popp, and J. L Cirac, Phys. Rev. Lett. 92, 027901 (2004). 

[11] A. Shimizu and T. Miyadera, Phys. Rev. Lett. 89, 270403 (2002). In this reference, states with p = 2 are called anomalously 
fluctuating states. 

[12] T. Morimae, A. Sugita and A. Shimizu, Phys. Rev. A 71, 032317 (2005). 

[13] A. Sugita and A. Shimizu, J. Phys. Soc. Jpn. 74 (2005) 1883. 

[14] S. Parker and M. B. Plenio, J. Mod. Optics 49, 1325 (2002). 

[15] Y. Shimoni, D. Shapiro and O. Biham, Phys. Rev. A 69, 062303 (2004). 

[16] Y. Shimoni, D. Shapiro and O. Biham, Phys. Rev. A 72, 062308 (2005). 

[17] R. Orus and J. Latorre, Phys. Rev. A 69, 052308 (2004). 

[18] R. Jozsa and N. Linden, Proc. R. Soc. London, Ser. A 459, 2011 (2003) 

[19] G. Vidal, Phys. Rev. Lett. 91, 147902 (2003). 

[20] A. Shimizu, talk presented at The 4th Symposium on Quantum Effects and Related Physical Phenomena (December 20-21, 
2000, Tokyo, Japan); A. Shimizu and T. Miyadera, in Proceedings of the 56th annual meeting of the Physical Society of 
Japan (Physical Society of Japan, 2001), paper no. 28pYN-6. 

[21] A. Ukena and A. Shimizu, Phys. Rev. A 69, 022301 (2004). 

[22] A. Shimizu and T. Morimae, Phys. Rev. Lett. 95, 090401 (2005). 

[23] T. Morimae and A. Shimizu, Phys. Rev. A 74, 052111 (2006). 

[24] For homogeneous states, it is obvious that A{L) should be homogeneously extended with increasing L. For example, if 

A{L) = 'Y^^=i(^z{l) for some L, then AilJ) for another L' (7^ L) should be X]^i<5"z(0- For this reason, the condition 
||a(Z)|| = 1 was not imposed in Refs. [HI, Fl^, [23I] . On the other hand, we do not require here that a{l) is independent of L, 
because the independence is not well-defined when treating general families of states. 

[25] R. Raussendorf, D. E. Browne, and H. J. Briegel, Phys. Rev. A 68, 022312 (2003). 

[26] D. Gross, J. Eisert, N. Schuch, D. Perez-Garcia, Phys. Rev. A 76, 052315 (2007). 

[27] This expectation has been confirmed for one-way quantum computers by Y. Matsuzaki and A. Shimizu, unpublished. 

[28] Y. Matsuzaki and A. Shimizu, in Proceedings of the 61st annual meeting of the Physical Society of Japan (Physical Society 
of Japan, 2006), paper no. 28aSA-4. 

[29] One can also apply the theories of Refs. [HI. [2^ to n-point correlations, if n (> 3) is independent of the system size L and if 
the n points can be divided into two groups 1, 2 in such a way that the maximum distance li (i = 1, 2) between two points 
in each group is independent of L. In fact, let I = max^^i, and regard each £ neighboring sites as a new site. The original 
lattice composed of L sites then becomes a new lattice composed oi L/l new sites. Then, the original n-point correlations 
become two-point correlations of the new lattice, and the theories of Refs. [HI . [2^ are applicable to the new lattice. By 
interpreting the results in terms the original sites, one obtains the corresponding results for the n-point correlations. 

[30] See, e.g., R. Haag, Local Quantum Physics (Springer, Berlin, 1992). 

[31] D. Deutsch and R. Jozsa, Proc.Royal Soc. London A 439, 553 (1992). 

[32] P. W. Shor, Phys. Rev. A 52 R2493. 

[33] J. Preskill, Proc. R. Soc. London, Ser. A 454, 385 (1998), and references cited therein. 
[34] N. Lynch, J. Assoc. Comput. Mach. 22, 341 (1975). 

[35] P. Orponen and U. Schoning, Mathematical Foundations of Computer Science 1984 (Springer, Berlin, 1984), pp. 452-458. 
[36] S. Even, A. L. Selman and Y. Yacobi, J. Assoc. Comput. Mach. 32, 205 (1985). 
[37] R. V. Book and D-Z. Du, J. Assoc. Comput. Mach. 34, 718 (1987). 

[38] Strictly speaking, this picture is not necessarily correct because when quantum states are strongly inhomogeneous Lemax is 
not necessarily proportional to the number of pairs of sites whose correlation is of 0(1). Fortunatelly, however, this picture 
is good for most instances in our simulations. 

[39] A. Carlini et al., Phys. Rev. Lett. 96, 060503 (2006), and references cited therein. 

[40] A. Hosoya, private communication. 



